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ABSTRACT 

Using the deep multi-wavelength MUSYC, GOODS, and FIRES surveys we construct a stellar mass-limited 
sample of galaxies at 2 < z < 3. The sample comprises 294 galaxies with M > 10 11 M Q distributed over four 
independent fields with a total area of almost 400arcmin 2 . The mean number density of massive galaxies in 
this redshift range p(M > 10 n M Q ) = (2.2±0.6) x I0~ 4 hj Mpc" 3 . We present median values and 25 th and 75 th 
percentiles for the distributions of observed /?ab magnitudes, observed J- K s colors, and rest-frame ultra-violet 
continuum slopes, M/Ly ratios, and U — V colors. The galaxies show a large range in all these properties. The 
"median galaxy" is faint in the observer's optical (Rab = 25.9), red in the observed near-IR (J—K s = 2.48), has 
a rest-frame UV spectrum which is relatively flat in F\ {(3 = -0.4), and rest-frame optical colors resembling 
those of nearby spiral galaxies {U-V = 0.62). We determine which galaxies would be selected as Lyman break 
galaxies (LBGs) or Distant Red Galaxies (DRGs, having J-K s > 2.3) in this mass-limited sample. By number 
DRGs make up 69 % of the sample and LBGs 20 %, with a small amount of overlap. By mass DRGs make 
up 77 % and LBGs 17 %. Neither technique provides a representative sample of massive galaxies at 2 < z < 3 
as they only sample the extremes of the population. As we show here, multi-wavelength surveys with high 
quality photometry are essential for an unbiased census of massive galaxies in the early Universe. The main 
uncertainty in this analysis is our reliance on photometric redshifts; confirmation of the results presented here 
requires extensive near-infrared spectroscopy of optically-faint samples. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The properties of massive galaxies at high redshift place 
important constraints on galaxy formation models (see, e.g., 
Kauffmann & Chariot 1998; Nagamine et al. 2005). The 
"standard" and most successful method for finding distant 
galaxies is the Lyman dropout technique, which relies on the 
strong break in the rest-frame ultra-violet (UV) spectra of 
high redshift galaxies blueward of the Lyman limit (Steidel 
et al. 1996, 1999). However, it is not yet clear whether these 
galaxies are representative of the high redshift galaxy pop- 
ulation, in particular at the high-mass end. As the Lyman 
break selection requires that galaxies are very bright in the 
rest-frame UV it may miss objects that are heavily obscured 
by dust or whose light is dominated by evolved stellar popu- 
lations. 

Advances in instrumentation have made it possible to se- 
lect galaxies in complementary ways, and recent studies have 
demonstrated that the universe at z > 2 is much more diverse 
than had been realized. Among recently identified "new" 
galaxy populations are submm galaxies (e.g., Smail et al. 
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2004), distant red galaxies (DRGs) selected by the criterion 
J-K s > 2.3 (Franx et al. 2003; van Dokkum et al. 2003), 
"IRAC Extremely Red Objects" (IEROs; Yan et al. 2004), 
and "BzK" objects (Daddi et al. 2004). 

The current situation is somewhat confusing, as the relative 
contributions of the various newly identified galaxy popula- 
tions to the stellar mass budget and the cosmic star formation 
rate are still unclear. Furthermore, as emphasized by, e.g., 
Adelberger et al. (2005) and Reddy et al. (2005), there can 
be considerable overlap between selection techniques, which 
makes it difficult to interpret the plethora of windows on the 
high redshift Universe. 

Ideally samples of high redshift galaxies are selected not 
by color or luminosity but by stellar mass. Whereas lumi- 
nosities and colors can vary dramatically due to star bursts 
and the presence of dust the mass evolution of galaxies is 
probably gradual. Also, galaxy formation models can pre- 
dict masses with somewhat higher confidence than luminosi- 
ties and colors. Stellar masses of distant galaxies are usually 
determined by fitting stellar population synthesis models to 
broad band photometry. Although there are significant sys- 
tematic uncertainties, the stellar mass of a galaxy is usually 
better constrained than the instantaneous star formation rate, 
age, or dust content (e.g., Shapley et al. 2001; Papovich et al. 
2001 ; van Dokkum et al. 2004; Forster Schreiber et al. 2004). 

In this Letter we explore the properties of a stellar-mass 
limited sample of galaxies. The main purpose is to measure 
"basic" aspects of massive galaxies to compare with simula- 
tions of galaxy formation: their density and colors. A sec- 
ondary goal is to quantify selection biases introduced by two 
of the most widely used techniques for identifying distant 
galaxies: the Lyman break technique and the J-K color se- 
lection of Franx et al. (2003). We assume Q m = 0.3, = 0.7, 
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FIG. 1 . — Relation between stellar mass and observed total K s magnitude 
for galaxies at 2 < Zphot < 3. The solid line shows our selection limit of M > 
10" Mq. The broken line shows the approximate photometric limit of the 
MUSYC survey, which accounts for the majority of galaxies in our sample. 
There are very few galaxies in FIRES and GOODS with M > 10" M e and 
K s > 21.3, and we estimate that the completeness of the full sample of M > 
10 11 M e galaxies is M 95 %. 

and Ho = 70 km s" 1 Mpc" 1 . All magnitudes are on the Vega 
system, unless identified as "AB". 

2. DEFINITION OF THE SAMPLE 

The sample is drawn from three deep multi-wavelength 
surveys, all having high-quality optical - near-IR photom- 
etry: the "ultra-deep" Faint InfraRed Extagalactic Survey 
(FIRES; Franx et al. 2003), the Great Observatories Origins 
Deep Survey (GOODS; Giavalisco et al. 2004) Chandra Deep 
Field South (CDFS), and the new Multi-wavelength Survey 
by Yale-Chile (MUSYC; Gawiser et al. 2006; Quadri et al. in 
prep). 

Photometric catalogs were created for all fields in the same 
way, following the procedures of Labbe et al. (2003). Pho- 
tometric redshifts were determined using the method of Rud- 
nicket al. (2001, 2003). Comparing the photometric redshifts 
with 696 spectroscopic redshifts (63 at z > 1 .5) gives a scatter 
in Az/(1 +z) of a = 0.06. Restricting the analysis to galaxies 
at z > 1.5 in the MUSYC fields gives a = 0.12, correspond- 
ing to Az w 0.4 at z = 2.5. Approximately 5 % of galaxies 
in this sample are "catastrophic" outliers. A full discussion 
of the quality of the photometric redshifts is given elsewhere 
(G. Rudnick et al., in prep). To determine masses, stellar 
population synthesis models were fit to the photometry using 
standard techniques (see, e.g., Shapley et al. 2001; Forster 
Schreiber et al. 2004). Bruzual & Chariot (2003) models 
were used, with solar metallicity (Shapley et al. 2004; van 
Dokkum et al. 2004) and a Salpeter (1955) initial mass func- 
tion from 0.1 Mq to 100 Mq. Star formation histories were 
parameterized by a declining star formation rate with char- 
acteristic timescale r = 0.3 Gyr (see Forster Schreiber et al. 
2004). The Calzetti (1997) reddening law was used, with 
extinction ranging from Ay = to Ay = 3. We note that the 
derived masses are probably not significantly affected by the 
presence of AGN (e.g., Rubin et al. 2004, Reddy et al. 2005), 
as their contributions to the broad band fluxes are probably 
small (Forster Schreiber et al. 2004; Webb et al. 2006). 



Fig-fflshows the relation between stellar mass and observed 
total K s magnitude for galaxies with photometric redshifts 

2 < z < 3. There is a clear relation, with a scatter of a fac- 
tor of ~ 10. We selected all 294 galaxies with 2 < z p h t < 3 
and stellar masses M > 10 11 M Q in the three surveys. The re- 
liability of this procedure was assessed in the following ways. 
First, we compared the masses derived from U-K s photome- 
try to masses derived from U-K s plus Spitzer/IRAC photom- 
etry in the HDFS (Labbe et al. 2005). Although the masses of 
individual galaxies can vary by ~ 30 % the systematic differ- 
ence is < 10 %. Next, we determined what fraction of massive 
galaxies are fainter than K s = 21.3, the approximate limit of 
the MUSYC survey. Only 5 % of galaxies with M > 10 11 M© 
in the deep FIRES and GOODS fields have K s > 21.3. Ex- 
tremely obscured massive galaxies could be missed even in 
the deep FIRES and GOODS data, but the fact that ~ 90 % 
of submm-selected galaxies at z ~ 2.2 have K <2l (Smail et 
al. 2004) implies that such objects are very rare. We conclude 
that our mass-limited sample of 294 galaxies at 2 < z < 3 is 
~ 95 % complete. 

3. DENSITY 

The FIRES, GOODS, and MUSYC surveys cover four in- 
dependent fields: FIRES MS 1054-03 (23 arcmin 2 ), GOODS 
CDFS 11 (69 arcmin 2 ), MUSYC SDSS 1030 (103 arcmin 2 ) 
and MUSYC HDFS (188 arcmin 2 ). The total area is 
383 arcmin 2 , of which 76 % is contributed by MUSYC. The 
average surface density of M > 10 11 Mq galaxies with 2 < z < 

3 is 0.71 arcmin" 2 , but there are large field-to-field variations. 
The density in the CDFS field is only 0.42 arcmin" 2 , 60 % of 
the mean and a factor of three lower than that of the highest 
density field, SDSS 1030. This large variation is indicative of 
strong clustering and implies that densities inferred from in- 
dividual ~ 100 arcmin 2 fields should be treated with caution. 

After a 5 % correction for incompleteness the mean space 
density p(M > 10 n M Q ) = (2.2±0.6) x 10" 4 Mpc" 3 . The un- 
certainty includes the effects of field-to-field variations, but 
does not include possible effects caused by systematic errors 
in the photometric redshifts (see, e.g., Shapley et al. 2005). 
We note that this density is a factor of ^ 5 lower than that of 
z fa 3 f/-dropout galaxies to R A b = 25.5 (Steidel et al. 1999), 
which typically have much lower masses. 

4. PROPERTIES OF MASSIVE GALAXIES AT 2 < Z < 3 

We use our mass-limited sample of 294 galaxies to deter- 
mine the median and dispersion in observed and rest-frame 
properties of the galaxies. Table 1 gives the median and 
25/75 -percentiles of the distributions of observed Rab mag- 
nitude and J—K s color; rest-frame U — V color and M/Ly ra- 
tio; and rest-frame UV slope, parameterized by F\ oc A' 3 . The 
rest-frame V magnitudes and U — V colors were determined 
from the observed magnitudes following similar procedures 
as those outlined in van Dokkum & Franx (1996). Rest-UV 
slopes (3 were determined from the best fitting spectral en- 
ergy distributions (SEDs), following the Calzetti, Kinney, & 
Storchi-Bergmann (1994) method of fitting to the ten rest-UV 
bins defined by those authors. 

As can be inferred from Table 1 the galaxies span a large 
range in all these properties. The "median galaxy" is red 
and faint in the observer's optical, with (Rab) = 25.9. We 
show the full distribution of the rest-frame U — V colors in 
Fig.|2ja). The bluest galaxies have U—V < -0.1 and are bluer 

1 ' Area with JHK S coverage. 
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FIG. 2. — Distribution of rest-frame U — V colors (a) and rest-frame UV 
slope P (b) for galaxies with M > 10" Mq and 2 < Zph ot < 3. The galaxies 
show a wide range in optical- and rest-UV colors. Blue histograms indicate 
galaxies with the colors and luminosities of LBGs; red histograms indicate 
DRGs. 
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FIG. 3. — Spectral energy distributions of three MUSYC galaxies with 
different rest-UV slope j3. Overplotted are the best fitting Bruzual & Chariot 
(2003) models. From top to bottom the galaxies have /3 = —1.3, —0.2, and 
0.9 respectively. The blue rest-UV SED of the top object is typical for LBGs. 
The slope of the middle galaxy is close to the median value of the full sample. 



than nearby irregular galaxies, whereas the colors of the red- 
dest galaxies are similar to those of nearby ellipticals (see, 
e.g., Fukugita, Shimasaku, & Ichikawa 1995). The median 
(U -V) = 0.6, which is similar to nearby spiral galaxies but 
also to nearby dust-enshrouded star burst galaxies (e.g., Ar- 
mus, Heckman, & Miley 1989). 
The distribution of (3 is shown in Fig.|2jb). Remarkably, the 



distribution is rather fiat and has no well-defined peak, in con- 
trast to previous studies of optically-selected samples (Adel- 
berger & Steidel 2000). The median (/?) = -0.39, indicating 
a relatively fiat spectrum in Fa. A potential worry is that in- 
dividual values of j3 are uncertain, as many galaxies are very 
faint in the observer's optical. We tested the robustness of 
the derived distribution of (3 by summing the observed optical 
fluxes of the galaxies in the lower and upper 25 % quartiles, 
weighting by the inverse of the total optical flux. The power- 
law slopes of these summed SEDs are in very good agreement 
with the median /3s that we determined from the SED fits. 

The large range of properties of massive galaxies at 2 < z < 
3 is illustrated in Fig. |3| which shows the full UBVRIzJHK, 
SEDs of three galaxies from the MUSYC survey with differ- 
ent values of (3. The top galaxy has a very blue SED similar to 
those of UV-selected samples (see, e.g., Shapley et al. 2001), 
the middle object has an SED which resembles that of nearby 
spiral galaxies, and the bottom galaxy has a very red SED in- 
dicating strong extinction. 



TABLE 1 

Observed and Rest-frame Properties 



Quantity 


25% 


Median 


75% 


iftoi.AB (obs) 


25.1 


25.9 


26.7 


J-K s (obs) 


2.22 


2.48 


2.85 


U-V (rest) 


0.41 


0.62 


0.80 


log M/Ly (rest) 


0.05 


0.20 


0.30 


(3 (rest) 


-1.13 


-0.39 


0.31 



5. DISCUSSION 

The main result of our analysis is that massive galaxies at 
Z ~ 2.5 span a large range in rest-frame UV slopes, rest-frame 
optical colors, and rest-frame M/Ly ratios, indicating signifi- 
cant variation in dust content, star formation histories, or both. 
This result is not surprising in the light of the recent discov- 
eries of DRGs, IEROs, and other populations. Here we have 
quantified the median colors and their range for a uniformly 
selected, large, mass-limited sample. 

The large variation in the rest-frame color distributions of 
our mass-limited sample implies that "standard" color selec- 
tion techniques produce biased samples. We consider two of 
the two most widely used selection techniques in this redshift 
range: the Lyman break technique of Steidel and collabora- 
tors and the J—K s > 2.3 Distant Red Galaxy (DRG) selection 
(Franx et al. 2003; van Dokkum et al. 2003). Lyman break 
galaxies are identified in the following way. From the best- 
fitting Bruzual & Chariot (2003) SEDs (which include absorp- 
tion due to the Ly a forest) we calculated synthetic colors in 
Steidel's U n G7Z system. To qualify as an LBG an object has 
to have Rab < 25.5 and synthetic U„G1Z colors which place it 
in the Lyman break, BX, or BM selection region (see Steidel 
et al. 2003, 2004). Combined these criteria provide a contin- 
uous selection of galaxies over the redshift range considered 
here. 12 Figure @] illustrates the LBG and DRG selection tech- 
niques, as applied to our sample. DRGs with J-K s > 2.3 
are indicated by red symbols and LBGs by blue symbols. 
The DRG limit and the "standard" photometric LBG limit of 
Rab = 25.5 are also indicated. 13 

12 A LBG in this definition is therefore an object which has Rab < 25.5, 
2 < Zphot < 3, and is either a classical "U-dropout" or a BX/BM object. 

13 We note that not all galaxies with J-K s > 2.3 have redshifts in the range 
2 < z < 3: to K s = 21 we find that ~ 50 % are in this redshift range with the 
rest about equally split between z < 2 and z > 3 galaxies. 
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FIG. 4. — Correlation between observed J —K s color and Sab magnitude. 
The majority of the galaxies are faint in R and red in J — K s . Red symbols 
denote DRGs with J —K s > 2.3; blue points denote LBGs with Rab < 25.5. 
Red symbols with blue circles fall in both categories. Massive LBGs have 
blue near-IR colors and are bright in the observed R band. The inset shows 
the optical color distribution of galaxies with R < 25.5. Only ~ 50% of 
optically-bright massive galaxies have the colors of LBGs. 

By number, DRGs make up 69 % of the sample and LBGs 
20 %. The DRG and LBG samples do not show much overlap: 
only 7 % of objects fall in both categories. By rest-frame V- 
band luminosity DRGs contribute 64 % and LBGs 32 %. By 
mass DRGs contribute 77 % and LBGs 17 %. Together, the 
LBG and DRG techniques identify 82 % of massive galaxies 
by number and 84 % by mass. Most of the remaining galaxies 
are optically faint, slightly bluer than the J-K s = 2.3 limit, 
and have redshifts z < 2.5. Approximately 85 % of them fall 
in the "BzK" selection region (Daddi et al. 2004), which is 
optimized for galaxies at 1.4 < z < 2.5. We note that the rel- 
atively small fraction of LBGs in the sample is not solely due 
to the imposed Rab <25.5 limit. As shown in the inset of Fig. 



|4|only — 50 % of galaxies with R AB < 25.5 have the rest-UV 
colors of LBGs, and this fraction decreases going to fainter R 
magnitudes: when no R limit is imposed we find that ~ 1/3 of 
the galaxies have the colors of LBGs. The underlying reason 
is the broad distribution of /?. 

It is clear from Fig.|4]that the LBG selection produces very 
different samples of massive high redshift galaxies than the 
DRG selection. Both samples are biased: LBGs are too blue 
and DRGs are too red when compared to the median values of 
the full sample. This bias is shown explicitly by the blue and 
red histograms in Fig.|2] The Lyman break criteria were de- 
signed to find star forming galaxies, but Shapley et al. (2004) 
and Adelberger et al. (2005) have argued that they can also 
be used to find massive galaxies at high redshift. This is ob- 
viously the case, but we find that the colors and M/Ly ratios 
of massive LBGs are not representative for the full sample of 
massive galaxies. Surveys over the full set of optical/near-IR 
passbands from U through K are essential to obtain represen- 
tative samples of massive galaxies. 

The main uncertainty in this analysis is the reliance on pho- 
tometric redshifts. We estimated the effect of this uncertainty 
by randomly perturbing the redshifts using a Gaussian distri- 
bution with dispersion Az/(1 +z) = 0. 12, and repeating the se- 
lection and analysis. Despite significant migration of galaxies 
in and out of the 2 < z < 3 redshift range the values in Table 
1 change by only ~ 10%. We note, however, that there may 
be subtle systematic biases which can have significant effects, 
in particular on the derived masses (see, e.g., Shapley et al. 
2005). Comprehensive tests of the techniques employed here 
and in other studies of infra-red selected samples (e.g., Dick- 
inson et al. 2003, Rudnick et al. 2003) are urgently needed, 
and will become feasible with the introduction of multi-object 
near-IR spectrographs on 8m class telescopes. 
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